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Abstract 

Techniques based upon the string organisation of amplitudes may be used to simplify 
field theory calculations. We apply these techniques to perturbative gravity and calculate 
all one-loop amplitudes for four-graviton scattering with arbitrary internal particle content. 
Decomposing the amplitudes into contributions arising from supersymmetric multiplets 
greatly simplifies these calculations. We also discuss how unitarity may be used to constrain 
the amplitudes. 



Address after Sept. 1, 1994: University College of Swansea, UK. 



1. Introduction 



Calculations in perturbative gravity are well known to be prohibitively difficult using 
conventional Feynman diagram techniques [1]. In typical gauges, the graviton vertices 
contain considerably more terms than in gauge theories and two powers of loop momentum 
rather than one. These two features help to make perturbative gravity calculation an 
algebraic nightmare. 

Recently, the alternate organisation of amplitudes offered by string theory has been 
used by Bern and Kosower to construct rules for computations in gauge theories [2,3]. 
These rules allow a considerable algebraic simplification compared to normal Feynman 
diagram techniques as evidenced by the first calculation of the five-gluon one loop am- 
plitudes [4]. Since string theory also contains gravity one can apply these techniques to 
obtain rules for calculations in perturbative gravity. In ref. [5], the string-based technique 
for perturbative gravity was outlined and a sample calculation of the ylone-ioop|-_^ _|_^ _|_^ 
four graviton helicity amplitude was performed. This is the simplest of the four graviton 
amplitudes being finite and without cuts and has been previously calculated by first cal- 
culating the contributions to the four graviton amplitude from real scalars in the loop and 
then using the supersymmetry Ward identities [6] . 

In this paper we present a detailed description of the rules for one-loop n-graviton 
amplitudes and use these to calculate the four graviton amplitude for all helicity configu- 
rations with arbitrary particle content in the loop. As with the QCD method, these rules 
arise by looking at the infinite tension limit [7,8] of a string theory. However, these rules 
can be used with no knowledge of string theory. The rules are in many ways a "double" 
copy of those for QCD, 

(Gravity) ~ (Yang-Mills)^ . (1.1) 

which reflects the fact that in string theory, a closed string may be regarded as the product 
of two open strings [9] 

(Closed String) ~ (Open String)^ . (1.2) 

This equivalence is largely true at the level of the integrands of diagrams. (There is, 
however, a small amount of interference between the two "halves" of the string which is 
related to the zero mode integral in string theory.) Using tree- level relationships which 
embody (1.1) [10] string theory has been used to calculate tree-level graviton amplitudes 
previously in ref. [11]. 

As inspired by the string-based method, we implement a supersymmetric decomposi- 
tion of the amplitudes similar to that recently used for gauge theories [12,13]. Instead of 
calculating the contributions to the loop amplitudes from individual particles circulating 
(which may be the graviton, gravitino, vector, Weyl fermion or scalar in a gravitational the- 
ory) we calculate the contributions from various supersymmetric multiplets plus the scalar 
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contribution. Corresponding to the five particle types we must calculate four supersym- 
metric contributions (which we choose to be the N = l, N = 4:, N = 6 and N = 8 matter 
multiplets) plus the scalar. The contribution from any individual particle type is just a 
linear combination of these. This decomposition enables us to exploit the simplifications 
found in supersymmetry calculations. In supersymmetric theories there are cancellations 
between the fermions and bosons. If one uses a suitable formalism, these cancellations are 
manifest diagram by diagram. This proves to be an enormous simplification. Examples of 
such beautiful formalisms are 1) the string based rules and 2) a superspace formalism using 
the background field method [14,15,16]. In a general gauge even in a superspace formalism 
the cancellations do not occur diagram by diagram. The relationship between string based 
rules for gauge theories and conventional field theory is particularly close when the field 
theory is organised using the background field method [17,18]. 

For the four-point function, without cancellations, the Feynman parameter integral 
generically has eight powers of Feynman parameters (or equivalently eight powers of loop 
momentum). For the N = 1 matter multiplet (containing a scalar and a fermion) there 
is a cancellation of the two leading powers of Feynman parameters which simplifies the 
calculation considerably. With increasing N more cancellations occur until for the maximal 
case, N = 8, all eight powers cancel and one is left with a trivial sum of scalar box integrals. 
The special cases of N = 8 supergravity and N = 4 super- Yang-Mills four point functions 
were obtained previously in ref. [19] also be examining the infinite string tension limit. 

Another advantage of the supersymmetric decomposition is that the supersymmetric 
amplitudes, again due to the cancellations in loop momentum, can be strongly constrained 
by unitarity via the Cutkosky rules [20]. In ref. [21,22], situations in a gauge theory 
where unitarity completely determines the one-loop amplitudes are given. For the four- 
point gravity calculation the N = 8 and N = 6 contributions can be determined completely 
using the Cutkosky rules and are in agreement with the string-based results. The Cutkosky 
rules are also used to check the cuts in the remaining amplitudes. As a further example of 
the uses of the Cutkosky rules we calculate the logarithmic part of the two-loop amplitude 

^2-loop(^ + ++). 

As is well known, pure gravity is renormalisable at one-loop [23] whereas gravity 
coupled to matter is not [24] . However, the ultra-violet infinities do not arise in ampli- 
tudes containing only external gravitons but appear in amplitudes with external matter 
[25]. Our amplitudes are all ultra-violet finite in agreement with the formal arguments. 
The amplitude with gravitons circulating in the loop contain the infra-red singularities as 
expected. 

2. Rules for one-loop gravity 

In ref. [2,3] rules were introduced for the calculation of gauge theory amplitudes. These 
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were obtained by taking the infinite string tension limit of string theory amphtudes and 
can be used instead of Feynman diagram techniques. Although derived from string theory, 
they can be used without explicit knowledge of string theory. The various contributions 
to the amplitude are associated with (p^ diagrams. Typically, the string organisation lead 
to a more compact integrand for these diagrams than that arising in conventional field 
theory. The string-based technique has been used to perform significant calculations such 
as the five-gluon one- loop contributions [4]. 

There are two slightly different but equivalent formulations of the rules; One is ob- 
tained by taking the infinite tension limit of a superstring [2,3] whereas the other is obtained 
by taking the infinite tension limit of a bosonic string [26,12]. The bosonic form of the 
rules is more compact but the appropriate rules for fermions circulating in the loop must 
be inferred from the superstring case. The string-based rules for one-loop gravity which 
we present here were outlined in ref. [5]. These are based upon the bosonic formulation. 
The rules for gravity have many similarities to those for gauge theories so we will be brief 
in presenting them. 

The initial step in the rules is to draw all (f>^ diagrams, excluding tadpoles. There is 
also no need to include diagrams with a loop isolated on an external leg since these vanish 
when dimensional regularisation is used. The external legs of these diagrams should be 
labeled, with diagrams containing all orderings included. The inner lines of trees attached 
to the loop are labeled according to the rule that as one moves form the outer lines to 
the inner ones, one labels the inner line with the same label as the most clockwise of the 
two outer lines attached to it. The contribution from each labeled n-point 0^-like diagram 
with legs attached to the loop is 



Jo JO 

X 



ni—2+e 



where the ordering of the loop parameter integrals corresponds to the ordering of the 
n£ lines attached to the loop, Xij = Xi — Xj. The Xi^ are related to ordinary Feynman 
parameters by Xi^ = Yl%=i '^i' -^red is the "reduced kinematic factor", which the string- 
based rules efficiently yield in a compact form. The lines attached to the loop carry 
momenta Pi which will be off-shell if there is a tree attached to that line. The dimensional 
regularisation parameter 2e = 4 — D handles all ultra-violet and infra-red divergences. 
The amplitude is then given by summing over all diagrams. We also use the equivalent 
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Schwinger proper-time form of the amplitude 



V 



s = I 



. n£ — l n 



^ ' -^1=1 iKj"^^ Km 



^ -^red {•^h ) • • • ) ^in^ ' ^) 



(2.2) 



as discussed in ref.[17] 

In order to evaluate K^ed, one starts with the graviton kinematic expression 



/n n 
dxidxi JJ^ exp 



ki 



exp 



X exp 



{ki ■ £j kj ■ £i) Gj^ 6' 



^3 



kj 



exp 



multi — linear 

(2.3) 



where the 'multi-linear' indicates that only the terms linear in all Si and Si are included. 
The graviton polarization tensor is reconstructed by taking efe^ — > ef^. From a string 
theory perspective is the bosonic Green function on the string world sheet, Gb and 
Gb are derivatives of this Green function with respect to left- moving variables, and Gb 
and Gb are derivatives with respect to right- movers. (Since a closed string is periodic the 
variables described the string world sheet can split into "left-moving" and "right- moving".) 
The term H^^ is the derivative of the Green function with respect to one left mover and 
one right mover variable. The functions and H^^ are to taken as symmetric in the 

i and j indices while is antisymmetric. Although the above expression contains much 
information in string theory, when one takes the infinite string tension limit [2,3] it should 
merely be regarded as a function which contains all the information necessary to generate 
iiTj-ed for all graphs. The utility of the string based method partially lies in this compact 
representation (which is valid for arbitrary numbers of legs!). The existence of an overall 
function which reduces to the Feynman parameter polynomial for each diagram is one of 
the most useful features of the string based rules. 

The appropriate expression for gauge theories is obtained by setting e = in the 
above. The gravity expression is like a double copy of the gauge theory expression apart 
from the H^^ terms which mix the left and right movers. 

The first step in applying the rules is to remove all of the G^ and G^ by integrating 
the kinematic expression by parts with respect to the variables Xi and Xi where necessary. 
When manipulating this formula we take G^, G^, G^ and G^ to mean dxiG^B^ 9^.G^, 
dxiG^B ^iiG^B respectively. (After direct substitution of the values of the functions 
in the field theory limit these relations are almost but not quite true; however, for the 
purposes of manipulating eq. (2.3) this distinction is unimportant.) While carrying out 
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this process one must take into account the cross-terms where a left-mover derivative hits 
right-mover terms, and vice versa. This can be done by using the results 



dxk 



d ■■■■ 



dxk ^ 

For example, if the expression 



exp 



1=1 



i<j 



ki • kj G ^ 



'^B y^B ) 



is integrated by parts with respect to xi, the result is 

4 4 



/ Y\ dxidxi exp 

i=l i<j 



ki • kjG ^ 



A34A12/7=?34\2 

'^B ^B yy'B ) 

x({ki- k^G^i + ki ■ k^G^i + fci • A;4Gij4)^i3 ^ ^ 



Having carried out the integration by parts we now may carry out simple substitution 
rules for each diagram to obtain JC^ed- First the nr<j ®^P[^« ' term and the inte- 

grals over Xi and Xi are dropped from the kinematic expression. (Since the appropriate 
contributions have been included in the rules). After integration by parts, /C will be a sum 
of terms each of which has n Gb and n Gb- (An Hb is equivalent to one Gb and one 

Any diagram will be a loop with ng legs attached with possible non-trivial trees at- 
tached to the loop. The rules have two parts. Firstly tree rules are applied to /C. These 
produce a truncated /C which is a series of terms each with ni Gb and Gb- Secondly 
loop substitution rules are applied which give the Feynman parameter polynomial for the 
diagram. The tree rules are applied iteratively working from the outside of the attached 
trees towards the loop. For a two-point tree with outer legs labeled by i and j, one carries 
out the substitutions 



2ki ■ kj ^ 



(2.4) 



in remaining factors 



in each term. This should be applied at each tree vertex. 

Once the tree rules have been carried out for a diagram, one applies the loop rules. 
These depend on the particles circulating in the loop. They are essentially independent 
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applications of the Yang-Mills rules to the left- and right- mover parts, with an extra 
substitution for cross-terms Hb- 

One-loop amplitudes depend upon the particle circulating in the loop, and the sub- 
stitution rules are corresponding different for different particle types. For gauge theories 
there are three types of particles/rules, those for scalars S, fcrmions F and vectors V . For 
gravity, the rules are a double copy of the gauge theory rules and this choice of substitutions 
can be chosen differently for the two copies. That is, we can apply different substitution 
rules to the Gb (left- movers) and the Gb (right movers). The particle content circulating 
in the loop corresponding to these choices of the loop substitution rules is given in table 1. 



Substitution 


Particle Content 


2[^, S] 


complex scalar 


-2[S,F] 


Weyl Fermion 


2[S, V] 


Vector 


-A[V,F] 


gravitino and Weyl Fermion 


mv] 


graviton and complex scalar 


A[V,V]-2[S,S] 


graviton 


-4[F,F] + 2[5, F] 


gravitino 



Table 1: Applying the substitution rules shown corresponds to having the particle 
content shown circulating in the loop, [x, y\ denotes applying substitution rules x and 

y to Gb and Gb- 

F and V each produce two types of contribution. The first contribution is just the 
scalar S but the second is different in the two cases. The different contribution we refer to 
as the "cycle" contribution Gy and Gp- 

F=S + Gf 
V^S + Gv 

The common S contribution is obtained by making the substitutions 

sign.{xij) + 2xij) 

sigii{xij) + 2xij) (2.6) 

in the Schwinger parameterization (2.2). (Before taking the infinite tension limit a S- 
function exist in ff^ however as discussed in ref. [8] this 5-function does not contribute 
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(2.5) 



G% - -^i- 

H% — 
B 2T 



in the infinite string tension limit of physical amplitudes.) The cycle contribution comes 
from "cycles" oi Gb- A cycle is a sequence of Gs's 



B ^B •••^B 



The substitution rules for these cycles is different in the three cases. For the scalar they 
are vanishing. For Cy, the substitution rules are 



^B '^B ^ ^ 

G'pG'p ...G'^-^'-'GY' ^1/2 (m>2) 

where all the cycles must follow the ordering of the legs, and only one cycle at a time 
may contribute to any term. Once these substitutions have been made all remaining G^'s 
should be replaced as in eq. (2.6). For Cp the following substitution is made 

m 

GTGT ■ ..GT^'-GT' - -(-1/2)- n ^^9n{x,.^.^,) 

k=l 

In contrast to the V rules, all cycles contribute in the F case regardless of ordering. Also, 
all combinations of one or more cycles from each term contribute. Again, once these 
substitutions have been made all remaining G^'s should be replaced as in (2.6). 

For example if, for the four-point amplitude, we have a term in the Kinematic expres- 
sion 

K = {G]if{Giy (2.7) 

Then for the box diagram with ordering of legs 1234 the cycle contributions for the two 
cases are 



Cf : K 
Cv : K 



1/1 \21/1 \2 1 

-5(3(1+2x34)) --(-(1 + 2X,.)) +- 

/I \2 /I n2 

(-(1 + 2X34)) +(-(1 + 2X12)) 



1^ (2.8) 



there being no cycle contribution in the scalar case. 

This process gives an expression for Kred for each diagram for arbitrary particle con- 
tent in the loop. The integral in (2.1) can now be carried out. The contributions from each 
diagram are then summed over. Explicit simple examples of the applications of the string 
based rules for QCD are given in refs. [3,12] and for gravity in ref. [5] which the interested 
reader may wish to examine to see the simplicity of the string based method. The string 
based rules have advantages in producing compact expressions for the numerators in the 
Feynman parameter integrals. 
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3. Supersymmetric Decomposition 



A useful way to organise the n-graviton amplitudes is to use a "supersymmetric de- 
composition" . The loop amplitudes depends upon the state circulating in the loop. In a 
graviton scattering calculation this may be one of five states: a scalar, a Weyl fermion, a 
vector, a gravitino or the graviton itself. The string based rules can be used to calculate 
these contributions individually. However, it proves convenient to calculate the contri- 
butions from supersymmetric multiplets instead. Amplitudes for all choices of particles 
in the loop can be written as linear combinations of those for certain choices of super- 
symmetric multiplets and for a scalar. In particular we choose one multiplet from each 
of N=l,4,6,8 supersymmetric theories with particle content given in the following table. 
These multiplets are centered around the spin-0 complex scalar. 



N 


scalars 


spin- 1/2 


spin-1 


spin-3/2 


spin- 2 


N = 


1 










N = l 


1 


1 








iV = 4 


3 


4 


1 






N = 6 


10 


15 


6 


1 




N = 8 


35 


56 


28 


8 


1 



Table 2: Particle content of the supersymmetric multiplets we consider. Scalars are 
complex, and the fermions are Weyl. 

This decomposition is useful when evaluating amplitudes. In general, the integrations 
involved in amplitudes increases considerable with the degree of the polynomial in the 
numerator of the Feynman parameter integral. (Or equivalently with the degree of the 
loop momentum polynomial if performing momentum integrals.) If one uses the string- 
based rules cancellations due to supersymmetry occur within each diagram, reducing the 
complexity of computations. Similar cancellations occur if one uses a background field 
method within a superfield formalism [16]. (The relationship between background field 
methods and string based calculations is explored in [17].) 

Specifically, within the string based rules these simplifications can be seen as can- 
cellations between the common contributions within multiplets. In terms of Feynman 
parameters, for a general n-point integral the scalar term is a polynomial of degree n 
however the cycle contributions are polynomials of degree n — 2. For each particle type, 
there is a scalar contribution Ng [S, S] where Ng counts the degrees of freedom with fermions 
having negative weight. Hence for any supersymmetric multiplet the [S, S] term will can- 
cel and the Feynman parameter polynomial will be simplified. With increasing N there 
are increasing cancellations. Also for the combination Cy — '^Cp the two- and three-cycle 
contributions cancel leaving a polynomial of degree n — 4. This may be seen, for example. 
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by comparing the cycle contributions in eqn. (2.8). For the supergravity multiplets can- 
cellations can occur on both left and right movers. For example, for the N = 8 calculation 
we have 



^A''=8 ^graviton g^gravitino _j_ 2g^vector gg^fermion _j_ gQ^scalar (3 1) 

Inserting the rules from table 1 we find 

^jv=8 ^ 4[(0^^ (j^^ _ ^2[Cy, Cf] + 6^Cf, Cf] = ^Cy - 4Cp, Cy - 4Cf] (3.2) 

From this we see that for a n-point integral the Feynman parameter polynomial would 
be 2n — 8 at most. The cancellations for a given are shown in table 3. (We are 
using a regularisation scheme which preserves supersymmetry [27] which simplifies the 
decomposition. In other regularisation schemes the form is the decomposition is a little 
more complex.) 



N 


Contribution 


Degree 


N = 


2 [5, S] 


2n 


N = 1 


2[Cf,S] 


2?i - 2 


N = 4 


2[Cv-4Cf.S] 


2n-4 


iV - 6 


-4[Cv -4Cf,Cf] 


2n-6 


N = 8 


A[Cv -4Cf,Cv -4Cf] 


2n-8 



Table 3: The String rules appropriate for the multiplet are given and the degree of the Feynman 
parameter polynomial for an n-point loop integral. 

To reconstruct the amplitudes for specific particles in the loop we can use 

^[0] ^ ^iV=0 



^[1/2] = ^N=l _ ^[0] 

^[1] ^ ^iV=4 _ 4^iV=l ^ ^[0] 
^[3/2] ^ ^Ar=6 _ g^Ar=4 ^ g^iV=l _ ^[0] 

^[2] ^ ^Ar=8 _ 8^iV=6 ^ 20A^=4 _ IQ^N=1 ^ ^[0] 



(3.3) 



4. Four-graviton amplitudes 

We now present the one-loop 4-graviton results for all choices of helicity. We will give 
each result in the supersymmetric decomposition form as described in the previous section. 
We also quote the pure gravity results explicitly. 
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The first step is to insert spinor helicity simplifications into the kinematic expression 
(2.3). The spinor hehcity method for gravitons [11,28] is related to that for vectors [29] by 



e'^'^ = e'^s^, s = e e 

where are the graviton helicity polarizations and are the vector helicity polariza- 
tions defined by Xu, Zhang and Chang. We use the notation for spinor inner products 
{k-\k+) = (12) and {k+\k-) = [12]. The use of spinor helicity techniques has proved 
extremely useful in QCD calculation. All states are taken to be outgoing and may have 
plus or minus helicity. There is no concept of colour ordering which is found in QCD 
amplitudes. There are thus three independent helicity configurations for the four point 
amplitude, (+, +, +, +), (— , +, +, +) and (— , — , +, +), the others being obtained by con- 
jugation from these. 

For the (— , -|-, -|-, -|-) and (-I-, -|-, -|-, -|-) one-loop amplitudes all the supersymmetric com- 
ponents in the decomposition vanish due to supersymmetric Ward identities analogous to 
the situation in QCD [30]. The tree level graviton amplitudes vanish for these helicity 
configurations and hence the one-loop results are the leading order for these configurations 
and have a simple form rather analogous to a tree amplitude without logarithms or infini- 
ties. From the inverse decomposition (3.3) this amplitude for any particle content is just 
proportional to the scalar contribution. The scalar contributions to graviton scattering 
have been previously calculated in ref. [6] but not in a spinor helicity basis. The results 
from the string-based rules agree with these results and we have checked explicitly that 
the cycle contributions cancel, demonstrating the Ward identities. We find 

Ail- 2+ 3+ 4+) - N j!^(^lY( ^il! Yil+^t+fl 

A[l ,Z ,6 ,4 i-^V.^^^^^U; l[12](23)(34)[41]7 5760 ^ ^ 

A(l+ 2+ 3+ 4+)--N ^( '-^ yis^ + st+fl 

A[l ,Z ,6 ,4 )- (4^)2 I, (12) (2 3) (3 4) (4 1)^ 1920 ^ ^ ^ 



where 



Ns = Nb- Nf 



is the number of bosonic states in the loop minus the number of fermionic states and 
s = {ki + k2)'^, t = (ki + k^)'^ and u = (ki + ks)'^. So, for instance, since a graviton is made 
up of two helicity states the amplitudes for pure gravity are found by putting A^g = 2 in 
the above expressions. 

For the ^(1~, 2~, 3+, 4+) amplitude, none of the cycle terms vanish, so these must be 
included. We express their contributions using the supersymmetric decomposition given 
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in the previous section. The (complex) scalar amplitude is 



rni, , F(t-u)(t^ + 9ut^ + 46uH^ + 9uH + u'^)\n(-t -u) 

^ ' ' ' ^ 30 s7 

F (2^4 + 23 ut^ + 222 uH^ + 23uH + 2u^) FuH^{ln{-t/ - + tt^) 
^ 180^6 ^8 

(4.3) 

where F is 

m^{4nYrr ( st {12)^ ^ _ z,sW(47r)-rr tree/.- ^- o+ u 
16(47r)2 V (12) (2 3) (3 4) (4 1);' 4(47r)2 ^ ' ' ' 

^tree(^;L-,2-,3+,4+) is defined in eq. (5.3) and 

rr = . (4.5) 

^ r(l-2e) 

The amplitudes for the supersymmetric multiplets given in table 2 are 

F{t^ + Utu + u^) Ft^u^(\n\-t/-u)+7r^) 

24s4 ^ 2s6 

F{t - u) {t^ + 8tu + u'^) ln(-t/ - u) 
1275 



A^=^ = ^(^{t- u)s M-t/ -u)- tu{\n^{-t/ -u)+ Tx^) + s"^ 



A 



,S2 



^Ar=8 ^ 2F / ln(-^) ^ In(-t) ^ In(-s) 



(4.6) 



e \ st su tu 

^ ^^ nn{-t)ln{-s) ^ H~u)H-t) ^ H-s)ln{-u) 
\ st tu us 

We chose to express the amplitude in the (unphysical) regime where all momentum vari- 
ables s, t and u are negative. One can obtain expressions in the physical region by the 
substitution 

In(-s) ^ ln(|s|) -Z7r0(s) (4.7) 

etc. 0(s) is the Heavyside function where ©(a;) = 1, a; > and Q{x) = 0,x <0. 

The pure gravity amplitude can be found using the expression in eq. (3.3) which gives 
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the result 



\€\st SU tu J 



2 ln{-u) In(-s) ^ 2 In(-t) In(-M) ^ 2 In(-t) In(-s) 



t-u 

(t + 2 (2 1 + (2 + 2 t^w - + 2 tw^ + 2 u^) {ln^{-t/ - u) + tt^) 

{t - u) (341 ^4 + 1609 t^-u + 2566 t^u"^ + 1609 tu^ + 341 u^) ln{-t/ - u) 

^ 30s^ 

1922 + 9143 t^u + 14622 t^u'^ + 9143 + 1922 u"^ 
^ 180^6 

(4.8) 

These expressions have the correct symmetry expected in the amphtude under, for example, 
interchange of legs 3 and 4. 

Only the N = 8 amplitude has 1/e singularities. These are purely IR singularities the 
result being UV finite. This is as expected since the N = 8 multiplet is the only multiplet 
containing gravitons and the loop amplitudes with other particles circulating are expected 
to be IR finite. 



5. Unitarity constraints &c consistency checks 

Unitarity, in the form of the Cutkosky rules, is a strong constraint on amplitudes. 
When cancellations occur in the loop momentum integrals it becomes particularly restrict- 
ing and if enough cancellations occur unitarity may be enough to uniquely determine the 
amplitude. Again, the supersymmetric decomposition is useful in this context since cancel- 
lations occur within supersymmetry multiplets. In ref [22] it was proved that, for a gauge 
theory amplitude, if the n-point loop integral has at most n — 2 powers of loop momentum 
in the numerator the amplitude is uniquely determined from the cuts. Unfortunately, since 
gravity amplitudes generally have 2n powers of loop momentum this result is not as useful 
in determining amplitudes since for large n the power of the loop momentum polynomial 
grows to be larger than n — 2 even in the case of maximal cancellation N = 8. In this case 
we may apply this result, in principle, for 2n — 8 < n — 2 that is for n < 6 to completely 
determine the amplitude. Looking at table 3 we see that for the four point amplitude 
this result may be applied to the N = 8 and N = 6 multiplet contributions. (Technically 
the result in [22] applies when the cancellations occur in the loop momentum polynomial 
rather than in the Feynman parameter polynomial, however by examining the zero mode 
integral, for example, in [17] one can see that a loop momentum representation exists for 
the string based rules and the result may be use.) In this section we will demonstrate how 
to apply unitarity to obtain the cuts in the one-loop amplitudes. For the four point N = 8 
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and N = 6 contributions the results agree completely with the string based results and for 
the remaining amplitudes the amplitudes are consistent with the cuts. 

To calculate the cuts in all channels using the Cutkosky rules. Consider the regime 
where one of the momentum invariants is positive and the remainder are negative. This 
allow us to isolate the cuts in a single channel. We can then use the Cutkosky rules to obtain 
the cuts. To apply the cuts one needs explicit, preferable compact, expressions for tree 
amplitudes. In general, we consider the cut in the channel {kmi+kmi+i+- ■ ■+A;m2-i+^m2)^ 
for the loop amplitude An-i{l, 2, . . . , n), depicted in fig. 1 and given by 

dLIPS(-£i,£2) At-^(-£i,mi,...,m2,£2) At--(-£2,m2 + l,...,mi-l,£i). (5.1) 



Instead of evaluating the phase-space integrals instead evaluate the off-shell integral 

^1 cut 



/ ^''"(-^i'"'i'---'"'2,^2)^ At^^^(-£2,m2 + l,...,mi-l,4)- 



(5.2) 



whose cut is (5.1). This replacement is only valid in this channel. In evaluating this 
off-shell integral, we may substitute = £3 = in the numerator; any terms with or 
£2 in the numerator cancels a cut propagator leading to an integral without a cut in this 
channel. Evaluating these cuts requires the tree amplitudes for all possible intermediate 
states, preferably in a compact form. For the four point those tree amplitudes which have 
been calculated previously in refs. [6,31] and in a helicity basis by Berends, Giele and 
Kuijf in ref [11] are sufficient. For pure gravity, the tree amplitudes A^'^'^^{ — h -|— 1-) and 
^tree^_l_ -|- -|— |-) vanish but A^^^^{ 1-+) is non-zero and is given by 

-^'"''^'■^'■^"■^")4( (12)(23)(34)(41) )^><! («-^) 

We first note that the A^-'°°p(-, -F, +, +) and A^-^°°p(+, +, +, +) one loop amplitudes 
have no logarithms and hence no cuts. This is consistent with the fact that there are no 
non-vanishing pairs of tree amplitudes which could contribute to the Cutkosky rules. 

To calculate the cuts in A^~^°°^{1~ , 2~ , 3"*", 4+), first consider the cut in the s-channel, 
{ki + k2)'^ as shown in fig. 2a. 

J (iLIPS^*'"^^(l-,2-,4,£+) X yl^'^^^(£^,£2,3+,4+) (5.4) 



z 
2 



This is non-zero for the case where the intermediate (cut) states are gravitons, however 
when the intermediate states are otherwise the tree amplitudes are zero. This is because the 
graviton vertex does not flip helicity of the fermions [11,32] hence the amplitude with two 
gravitons and two fermions of the same helicity vanishes, A{g , g , ip'^ , ip'^) = 0. Similarly 
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the vector amplitudes and scalar amplitudes (with the concept of helicity being replaced 
by particle/antiparticle) vanish. Hence states other than gravitons do not contribute to 
eq. (5.4). For the supersymmetric decomposition this implies the s-channel cut will only 
be non-zero for the N = 8 contribution. Inserting the graviton tree amplitudes into (5.4) 
yields, 

32 J V (12) (2 £2) (^2 4) (41); (/ci-4) V(34)(44)(44)(43)y' {h ■ £2) 

(5.5) 

which we can rearrange, using the fact that £1 and £2 are onshell, 

32' [{12){34)J J ''''''^\h.£2){k,.£2){2£,f{£,lf{A£,f{£,3f ^ ^ ^ 

We can rearrange this, using 

{£i£2f (4 £2)' [12]' 



{£rlf{2£2f (41) (24) (41) [12] [12] (2^2) 

(4 4)' [12]' 



(41) (24) (44) [42] [14] (44) 

[12]^ 

(24) [42] [14] (41) 

[12]' 



4(/c2-4)(/ci-4) 

where we use {£1 1) [12] = - {£1 £2) [£2 2] etc. Similarly 

(4 4)' [3 4]' 



(44) (43)' 4(A;3-4)(/c4-4) 

Noting that 

\^ 91 \?, 41 = 

(2 3) (41) 



we obtain the form of the cut 

in'^ .tree/-,- o- o+ . + N„3.„. /" ,t ttoo (^2 " 4) (fe " 4) ^, 1 



;^tree(^-^2-,3+,4+)s^t« j dLIPS^ 



X 



(5.7) 



(5.8) 



[12] [3 4] = —^ (5.9) 



8 ^ ' ' ' ' y (fci-4)(fc4-4) i6(fc2-4)(fci-4)(fc3-4)(A;4-4) 

Q. f 1 



'.A^ree^-^- ,3+ ,A + )sHu j dLIPS 



16(A:i-4)(A:4-4)(A;i-4)(A;4-4) 

(5.10) 
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There is a useful identity, requiring the fact that both trees are on-shell. 
1 



-11 1-1 

+ Ti r-TT-, + T, ^-TT-, + 



s''\{k^-h){k4-h) {ki-£2){k4-h) {k^-h){.k4-h) {k^ ■ h){.k4 ■ h) 
16 . 



s^\{ki-£ink4+hr {ki+£2nk4+£ir {ki-hnk4-£2r {ki + £2^4 - £2^ 

(5.11) 

This uses momentum conservation extensively. Inserting the two propagators I /if and 
1/^2 and replacing J dLIPS by J d^i/{2n)^ as in eq. (5.2) the cut in eq. (5.10) can now 
be recognised as the cut of the sum of two scalar box integrals with orderings 1234 and 
2134 ( the four terms above only correspond to two independent boxes.) These boxes have 
coefficients, 

2^tree^^-^2-,3+,4+)stw (5.12) 



Next, consider cuts in the channel {ki + k^)'^. These in general are more complex and 
also depend upon the multiplet under consideration. As can be seen from fig. 2b in this 
case all particles contribute. To evaluate the cut one needs the four point amplitudes with 
two external gravitons and two scalars or fcrmions or vectors or gravitinos. These may be 
obtained from the four graviton amplitudes using an extended form of the supersymmetric 
ward identities [30] . From these, we obtain 

^(^-,0-,0+,^+) = ^l^A{g-,g-,g+,g+) 

(1 2) 

Aig-,A-,A+,g+) = !^A{g- , g' , g+ , g+) 

^12)^ (5.13) 

A{g-,A-,A+,g+) = ^^A{g-,g-,g+,g+) 

(1 2} 

Aig-,^-,^+,g+) = ^^^Aig-,g-,g+,g+) 

where (^r, ip, A, A, (p) are the members of the N — 8 multiplet. 

It is useful to count each states contribution to the cut relative to that for a scalar. 
All these states will contribute to the cut, 

j dLIPSA(4+,l-,(/>(£2)-,</>(^i) + )xA((/>(£i)-,(/>(£2)+,2-,3+)xpw=8 (5.14) 
which is explicitly 

f (l^i)' (2^i)' (1 4)^ (24)^ {ki ■ 4)(A:2 • ^2) ..... 



32 (23)^(41)^7 (3£i)^(4£i)^(£i^2) (^i • ^i)(^2 • 4) 
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The factor p for AT = 8 will be 



pN^8 = x^ - 8x^ + 28a;^ - 56x'^ + 70 - BGx'^ + 28x-^ - 8 * x"^ + x'^ (5.16) 

where 

, _ (U,) (2 4) 

The central term, 70, arises from the 35 complex scalars, the x^ and x~'^ from the Weyl 
fermions, the x^ and a;~^ from the vectors, the x^ and x~^ from the gravitinos and the x^ 
and x"^ from the gravitons. These relative weights compared to the scalar contribution 
are obtained from eq. (5.13). This simplifies to 

. -1^8 (^'-1)' ((l£2)(2£i)-(14)(2£2))^ 

PN=8 = [X — X ) = — 



X^ (l£2)(2£l)(l£l)(2£2))^ 

(12)'(£l£2)' 



(5.18) 



((1^2) (2£i)(l4)(2£2))^ 

using the identity (a 6) (cd) — (ac) (bd) + (ad) (cb). The cut then becomes 
i^%2 (12)^ /^TTpq (fei-^2)(fe2-^2) (4^2)^ 

32 (23)'(41)'y (A;i-£i)(A;2-£i)((l£2)(44)(£22)(3£i))2 ^ 

Similarly to before we have the identities 



{hhf [41]' {iii2f [2 3]' 



(£i 4)' (1 e^f 4(A;i • £2)(A;4 • h) ' (3£i)' (£2 2)' ■ h){k^ ■ h) 

and 



(5.20) 



[41]' [2 3]'= ,_f2l..2 (5-21) 



(12)' (34) 
which give the cut to be 

32 V (1 2)' (2 3)' (34)' (4 if) J (h ■ h) {k2 ■ h) 16(A;i • i2){k^ ■ h){k2 ■ l2){kz ■ 4) 

(5.22) 

which is 

This is just the analogue of the previous case. We can carry out the same factorization 
as before and find the two appropriate boxes with ordering 1234 and 1243. The cuts in 
the {k\ + kzf channel are obtained exactly as the t-channel cuts. Replacing, J dLIPS by 
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J d^£/{27r)'^ we can thus deduce that the cuts in all channels are described by the sum 
over the scalar boxes with coefficients in'^A^'^^^stu/A. The scalar box, for ordering of legs 
1234 is 

rri^nY 1 f 2 



h = 



(47r)2 st\e 
rr{4nY 1(4 2 



(47r)2 stie2 e 
This sum over boxes then evaluates to 



(ln{-s) + In(-t)) + 21n(-s) In(-t) - 7r^| 



(5.24) 



L^=«(l-, 2-, 3+, 4+) = in^!l^A'^^%l-,2-,3+, 4+) 



{4n) 



,'2sln(-s) + 2tln(-t) + 2Mln(-t() , . . . , . , , , . x, . x 

X ^ ^— ^ — ^— ^ + 2s ln{-t) In(-'u) + 2t In(-'u) In(-s) + 2u In(-s) In(-t) 



(5.25) 

Since we have produced an expression with the correct cuts which is written in terms of 
integral functions we can use the results of rcf. [21,22] to deduce that this expression is 
the entire amplitude. This is in agreement with the explicit calculation of the previous 
section. 

For the other cases, the, = 6,A^ = 4,N = 2 and = contributions, we may 
perform a similar calculation to the N — 8 case but with differing p. We will have 

PN=6 = [x -X ) - 



3 



((l£2)(2£i)(14)(2£2)) 

U .-1^4 (12)" (^1^2)" (roa) 



PN=2 = (X-X ) = 



((1^2) {2£i) {l£i) (2^2)) 
A similar calculation of the cuts may be performed. For the N = 6 calculation the cuts are 
enough to completely reconstruct the amplitude and we obtain an amplitude in agreement 
with the explicit string based-calculation. For the other case, the cuts would not be 
enough to completely specify the amplitude but provide strong consistency checks on the 
amplitudes. In general, the cuts are relatively simple to calculate when compact forms for 
the tree amplitudes exist and where the cuts specify the amplitude completely they would 
be the calculational method of choice. 

In the next section we will illustrate the use of the cuts to obtain the logarithmic parts 
of two-loop amplitude A{+, +, +, +). 

6. The Cuts in vl^-loop^^^ ^) 

In this section we will use the Cutkosky rules to calculate the cuts in the two-loop pure 
gravity amplitude A'^~^°°'^{+, +, +, +) . It is possible to do this because of the simple form 
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of the one-loop amplitude A^~^°°p {+,+,+,+). If we consider the two-loop cuts we must 
consider the cuts in the three-particle intermediate states as shown in fig. 3a. Fortunately, 
these cuts vanish because one of the two tree amplitudes must inevitably have a single 
negative helicity and this tree vanishes. Considering the remaining possibilities we find, 
for example in the s-channel that the configurations in figs. 3b and 3c may contribute. In 
this case we have a product of , — , +, +) and +, +, +) neither of which 

vanish. Since the one-loop amplitude does not contain logarithms or dilogarithms the 
evaluation of the cut is analogous to calculating a one-loop cut and we are able to do so. 
Explicitly the s-channel cut in fig. 3b is 

J dLIPS(-£i,^2) A'-'-°^{l+,2+, -£+,£+) A'^-%-£-,£-,3+,4+), (6.1) 

where ciLIPS(— £i, £2) denotes the Lorentz- invariant phase space measure. The explicit 
form of the tree and one-loop amplitudes is obtainable from eq. (4.2) and eq. (5.3). With 
the parameterisation shown 



^tree, ^- ^- 3+ 4+^ _ ^^V (^2 h)' V s{2k: 

[ t^,t,,6 ,4 j- ^ \{£^3){34){4£,)J 



3 ■ t2 



4 V(^2 3)(34)(4£i); {2k4-£2) 

(6.2) 

The evaluation of this cut employs many of the tricks already used to obtain checks 
of the cuts for the one-loop results. Using cq. (5.7) and (5.8) we can simplify the cut to 



1 ^5.2 r (2fc2-£2)f(s^ + s(2A;2-^2) + (2/C2-^2)^ 

^^LIPS .o,: ..... .^.o. .^ (6-3) 



1920 ((12) (2 3) (34) (41))2 J {2ki ■ £i){2k4 ■ £i){2k4 ■ £2) 

Using 

1/1 1x1/1 1 N / 



{2k4-£^){2k4-£2) s\2k4-£i) {2k4-£2)' s^{k4 + £i)'^ (/c3 + 4)2' 

this reduces to a sum of two box integrals. However since (2/c2 ■ £2) = — (2/ci ■ £1) the box 
integrands further reduce to the triangle integrals shown in fig. 4 with loop momentum 
polynomial 

(^s"" + s{2k2 ■ £2) + {2k2 ■ £2)^) (6.5) 
Evaluating this integral gives the result, 

1^^+^^ + ^(3«' + - 2s') + \{lu' + It' - 5.2)) X {-s)-^-^ (6.6) 
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The second triangle gives the same result and the contribution from fig. 3c is also the same. 
This equation has the correct cut in the s-channel. The t and u channel cuts are obtained 
from this by permuting s,t and u. We can thus obtain an expression for the logarithmic 
parts of A'^-^°°P{1+, 2+, 3+, 4+) 



y42-ioop(i+, 2+,3+, 4+) ~ Ai-i°°P(l+, 2+, 3+, 4+) 



X 



i-sy-' + i-ty-' + {-uy~' ^ i (sw^ + 3t^ - 2s^) {-sy- 



2e2 2 (s'^ + t^ + u^) e 

1 {3u^ + 3s^ - 2t^) i-ty-' 1 (3g^ + 3t^ - 2^^) {-uy-^ (6-7) 

^2 (§2 + ^2 + ^2) e +2 (s2 + ^2 + ^2) e 

+ polynomials 



where the polynomial pieces are not obtainable from the cuts. Associated with the log- 
arithms in the above expression are 1/e poles. The A{+ + ++) two-loop amplitude is 
an interesting object in perturbative gravity. Gravity is non-renormalisable at two loops 
however not all amplitudes contain obvious non-renormalisable UV infinities. As is well 

known [33], the two loop amplitude A{ does not contain such infinities and the 

infinities reside in the so-called "helicity-fiip" amplitudes A{+ + ++) and A{ — h -|-+). 
Calculation of the UV infinities for these amplitudes is a considerable undertaking and the 
explicit verification of the non-renormalisability of gravity at two-loops was an important 
but difficult calculation [34,35]. Although we have a suggested form for the infinities in this 
amplitude, obtained with an almost trivial calculation, we are not able to extract the UV 
infinity. Specifically the recognising of the UV and IR infinities is not possible by examin- 
ing the cut integrals. When the cut was reduced to a triangle integral with polynomial in 
eq. (6.5) one would normally extract the UV infinity from the which would give 

^ X 6^, (6.8) 

However this yields zero since the coefficient is k^k'(. Since the cut is only sensitive up to 
terms proportional to and £2 the loop momentum polynomial could have been replaced 

by 

((2A;i • £2)' + (2A;i • £2) (2/c2 ■ h) + {2k2 • £2)') (6.9) 

without affecting the cut. However this polynomial would give an UV infinity since the 
5-function would no longer vanish for the middle term. Thus although tempting, we are 
unable to deduce the coefficients of the UV infinities in the two-loop amplitude, at least 
without further information. 

For the configuration ^(1~, 2+, 3"'', 4+) the form of the one-loop amplitude also does 
not contain logarithms and one may evaluate the equivalent cut diagrams to those in fig. 3b 
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and 3c however in this case there is a diagram as depicted in fig. 5 which is non-zero and 
needs genuine two-loop integrals to evaluate. 

Two-loop amplitudes are formidable calculations in general. Progress towards a string 
based or string inspired method has been made both based upon the infinite string tension 
limit [36] and upon the world-line formalism [37] . In any formalism we expect unitarity to 
provide very useful checks upon the calculations as evidenced by the simplicity by which 
we obtained the cuts in ^2-ioop^_|_ _|_ _|__|_^ (albeit the simplest case possible). 

7. Conclusions 

Recently new rules for calculations of one-loop amplitudes have been constructed using 
string theory methods. These have been successfully used to do calculations in both gauge 
theory and gravity which have not been practical using conventional methods. In this 
paper we gave a detailed description of a set of rules for gravity derived from string theory. 
We used the rules to calculate one-loop amplitudes for four graviton scattering. These 
covered theories with arbitrary particles content. Our results are consistent both with 
previous calculations and formal arguments concerning divergences. Although perturbative 
quantum gravity is a non-renormalisable field theory the UV divergences do not appear in 
our calculations. 

In order to simplify the calculations we used a supersymmetric decomposition of the 
amplitudes also inspired by string theory. In a supersymmetric theory, provided one uses 
a suitable formalism, there are generally a large number of cancellations between different 
particle amplitudes. Careful choice of the supersymmetric multiplet amplitudes calculated 
enabled us to exploit these simplifications. Individual particle contributions could then be 
found from linear combinations of these supersymmetric amplitudes. The supersymmetric 
decomposition in the string based rules has proved extremely useful since it reduces the 
degree of the momentum loop polynomial diagram by diagram. Such cancellations do not 
occur even in a normal superfield formalism but are familiar if one uses a background field 
superfield formalism. The calculational advantages of calculating iS-matrix elements using 
background field methods are widespread especially in situations where many pure gauge 
vertices appear or in situations where cancellations are possible. 

We used unitarity constraints to check the amplitudes calculated. These constraints 
were found by use of the Cutkosky rules. For the N — 8 and = 6 multiplets unitarity 
determined the amplitudes completely. For the remaining cases those parts of the ampli- 
tudes containing cuts could be checked. Since the loop momentum polynomial grows as 
2n for an n-point amplitude rather than as n for gauge theories the Cutkosky rules are not 
as powerful a constraint in perturbative gravity as in gauge theories. 

We also showed how the Cutkosky method could be used to obtain the logarithmic 
parts of a two loop gravity amplitude by calculating the cuts for the two-loop +, +, +) 
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helicity amplitude. This amusing calculation is incomplete since there are potential poly- 
nomial terms but it does illustrate the potency of unitarity. 

In conclusion we have found considerable calculational benefits from calculating using 
the string-based methods [2,3] and derivatives thereof [21,22]. Many of the techniques 
motivated by string theory have proved useful in these gravity calculations and we expect 
these to have wider validity. 
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Figure Captions 

1 A generic cut in the amplitude is given by an integral over lorentz-invariant phase 
space (LIPS) of the product of two tree amplitudes. 

2 The helicity configurations for the s and t cuts in A^~^"°p[-\-., +, +, +) 

3 The three possible contributions to the s-channel cut for +, +, +). Con- 
tribution (a) is vanishing since there is no choice of helicity for the intermediate legs where 
both trees are non-vanishing. 

4 After manipulations, the s-channel cut in (+,+,+, +) reduces to a sum of 
simple triangle integrals with integrands quadratic in the loop momentum. 

5 This helicity configuration is non-vanishing and contributes to the cuts in ^2-Zoop^_^ _|_^_ 
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Fig.l 
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